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Based on the London approximation, we investigate numerically the stability of the elementary 
configurations of entanglement, the twisted-pair and the twisted-triplet, in the vortex-lattice and 
-liquid phases. We find that, except for the dilute limit, the twisted-pair is unstable and hence 
irrelevant in the discussion of entanglement. In the lattice phase the twisted-triplet constitutes a 
metastable, confined configuration of high energy. Loss of lattice symmetry upon melting leads to 
deconfinement and the twisted-triplet turns into a low-energy helical configuration. 
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The combination of the soft elastic moduli and the 
large temperatures attainable in the vortex system of the 
high-T c superconductors boosts the importance of fluc- 
tuations and leads to interesting phenomena such as vor- 
tex lattice melting and the appearance of vortex-liquid 
phases 0. In this context, topological excitations in the 
vortex system leading to entanglement of the flux lines 
play an important role, both with respect to statistical 
mechanics as well as dynamical properties of the vortex- 
solid and -liquid phases. In this letter, we present a de- 
tailed analysis of the stability and recombination prop- 
erties of the elementary entangled configurations, the 
twisted-pair and the twisted-triplet (see Fig. 1), for both 
the vortex-solid and for a model vortex-liquid phase. 

Topological excitations of the vortex-lattice in the form 
of edge- and screw-dislocations are long time known ob- 
jects Q|. Recently, interest has concentrated on more ex- 
otic configurations such as interstitials and vacancies || . 
The latter are relevant in the discussion of a novel super- 
solid phase in layered high-T c superconductors and 
can be viewed as bound pairs of oppositely "charged" 
edge dislocations. Similarly, vortex entanglement is con- 
sidered to be at the basis of yet another new vortex 
phase, namely the vortex- liquid Q, where the entangle- 
ment loops can be viewed as bound pairs of concentric 
screw dislocations of opposite sign. The role played by 
such basic loops of entanglement on the statistical me- 
chanics of vortex- liquids is the following 0] : If the equi- 
librium state of the vortex-liquid contains loops of all 
scales with a finite density, the liquid exhibits a normal 
(dissipative) response under application of a longitudinal 
current density j || B, and the entangled vortex- liquid is 
equivalent to the normal metallic phase ||. If, on the 
other hand, the vortex-liquid remains disentangled, lon- 
gitudinal superconducting response survives the melting 
transition and we obtain a new intermediate liquid phase 
distinct from the normal metallic one 0. 

Vortex entanglement is equally relevant for the dynam- 
ical properties of the vortex-liquid: the prohibitively long 



relaxation times via reptation |8[ make the barriers for 
vortex reswitching the limiting factor in the vortex dy- 
namics. The reswitching barriers then determine the in- 
ner viscosity of the liquid and hence its pinning, creep, 
and flow properties 

Previous work on vortex cutting [p|JlO|l has produced 
first estimates for the reswitching barrier of two isolated 
vortices; however, as we will show below the relevance 
of these results for the entanglement problem in vortex- 
solids and -liquids is rather limited. Wilkin and Moore 
have determined the excitation energy of the cross- 
ing configuration of a vortex pair in a vortex-lattice, as- 
suming this configuration to constitute a saddle for vor- 
tex reswitching. However, as we will show below, there 
is no metastable state for the twisted-pair configuration 
in a vortex-solid and hence there exists no saddle for 
reswitching (a similar result has recently been obtained 
by Dodgson and Moore |l2| within the framework of the 
lowest Landau level approximation valid close to H C2 ). 
A metastable twisted-pair configuration does exist in a 
model vortex-liquid, where the pressure of the surround- 
ing vortices acting on the pair is modelled by a circular 
potential. This situation has been investigated by Car- 
raro and Fisher [|l3) , who calculated the reswitching bar- 
rier for the limiting case of an infinitely extended twist 
using quite ingenious symmetry arguments. Below we 
will argue that in a realistic description of the vortex- 
fluid the pressure excerted on a vortex pair by its sur- 
roundings destroys the metastable twist in the same way 
as in a vortex lattice and hence entanglements involving 
three or more vortices have to be considered. 

Here, we study the stability and the reswitching barrier 
(connecting the metastable entangled state with the sta- 
ble rectilinear groundstate through a saddle) for the ele- 
mentary entangled configurations, the twisted-pair (TP) 
and the twisted-triplet (TT). The analysis is done for an 
isotropic superconductor (penetration depth A, coherence 
length £) within the London approximation; generaliza- 
tion of the results to the anisotropic situation involves 
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simple rescaling jjj , at least for B ^> H Cl . We start 
out with the vortex-lattice and show that no metastable 
twisted-pair configuration exists for fields B > H Cl . The 
TP state can be rendered metastable either by artificially 
enhancing the vortex core energy or by going over to 
the dilute limit B < H Cl , where the interaction between 
vortices becomes short range. Next, we consider the 
twisted-triplet in the vortex-solid and find a metastable 
state. The twist is restricted to a finite length along 
the field axis and we call this a "confined" excitation. 
The confinement is a consequence of the discrete lattice 
symmetry and leads to a high excitation energy when 
compared to the rectilinear groundstate. In comparison, 
the barrier stabilizing the TT state against reswitching 
is small and one concludes that the lattice phase shows 
only little entanglement. In the model vortex-liquid the 
situation remains essentially unchanged as regards the 
twisted-pair — the pressure of the neighboring vortices 
destroys the metastable state in the same manner as in 
the vortex-solid phase. For the twisted-triplet the as- 
sumption of a circular effective potential mimicking the 
liquid environment is realistic. The restoration of pla- 
nar rotational symmetry upon melting leads to decon- 
finement and the TT turns into a low-energy helical con- 
figuration stabilized by a high barrier against reswitch- 
ing. As a consequence, one expects the liquid to become 
entangled, however, due to the large reswitching barrier, 
a non-entangled system only slowly transforms into an 
entangled state. Whether a thermodynamic vortex-fluid 
phase becomes entangled immediately upon melting is 
a complicated statistical mechanics problem and we 
will not go into this discussion here. In the following, we 
give a brief description of the (numerical) technique on 
which our analysis is based and then present the results 
for the twisted-pair and twisted-triplet configurations in 
the vortex-lattice and the model vortex-liquid. 

We base our analysis on the London approximation 
valid for the important low and intermediate field range 
B < 0.2 H C2 . Choosing a set of n vortices (labelled by 
/i, v = 1, . . . n) involved in the excitation, the free energy 
functional takes the form (see Ref. JL|, r M = (R M ,z)) 



where r M and ry (r^ = r M — iy ) refer to separate points 
on the same line. Here, we treat the constant c describing 
the vortex core energy as a parameter; within the London 
model its physical value is c « 0.5. Taking all nearest 
neighbors into account (see Fig. 1), we choose n = 10 
and n = 12 for the TP and the TT, respectively. Twisted 
metastable configurations are obtained from a topologi- 
cally correct initial state with subsequent application of a 
conventional conjugate gradient method to minimize the 
energy. In case a metastable state exists, the reswitch- 
ing barrier is found by imposing a constraint dragging 
the configuration from the metastable minimum over the 
saddle towards the rectilinear groundstate. Within a con- 
strained configuration the vortex pair/triplet is forced to 
have a prescribed distance d in the z = symmetry 
plane. With this constraint imposed, the metastable min- 
imum is force-free, however, the saddle configuration in 
general remains forced and hence we obtain only an up- 
per estimate for the energy of the true (force-free) saddle 
configuration. Wc point out that neglecting the contri- 
bution of the core to the self-energy !F se i{ produces severe 
instabilities leading to unphysical fluctuations upon re- 
laxation. The origin of these short-range single-vortex 
fluctuations can be traced back to the dispersive nature 
of the line tension £i{k z ) = e hi(l/fc z £). The functional 
containing a finite core energy (c ~ 0.5) is not only 
physically correct but also stable with respect to such 
pathological fluctuations. 

Before turning to the specific discussion we briefly 
mention the natural scales in the problem: For a con- 
fined excitation the natural length scale along the field 
axis is the lattice constant a , whereas the scale for the 
excitation energy is e a . For the (uniaxially) anisotropic 
situation these scales change to ea Q for the length and 
ee a for the energy, where we have introduced the ef- 
fective mass ratio e 2 = m/M <C 1 M- The energy scale 
is conveniently expressed through the melting tempera- 
ture M, ee a « T m /2.7c^ ~ 6T m , where we have used 
a Lindemann number c L = 0.25 in the last equation. 

Vortex-lattice: For the interaction with the surround- 
ing vortices we choose the lattice potential 



with 



with R m denoting the equilibrium lattice sites and K 
is the zero-order modified Bessel function. For inter- 
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< B < H C2 we can use the limit 
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mediate fields H c 
X — > oo in V, leading to a simple and rapidly conver- 
gent series upon resummation in Fourier space. We then 
proceed along the lines described above and first search 
for a metastable twisted-pair configuration. For fields 
B > H Cl and using the correct core energy with c ~ 0.5 
no such metastable state exists; in fact, although such a 
state does exist for the isolated pair, it is squeezed away 
by the pressure of the surrounding vortices as the lattice 



2 



potential V is switched on. Thereby the twisted vor- 
tices tend to align antiparallel in the crossing region and 
the resulting attractive force between the segments leads 
to their collapse and mutual annihilation. This tendency 
of antiparallel alignment can be suppressed by artificially 
increasing the vortex core energy: for large enough c such 
that c > c c « 4.9 — ln(a /£) the metastable TP is recov- 
ered. Extrapolating this result (obtained for a /£, < 20) 
to smaller fields, a metastable TP configuration is pre- 
dicted for a /£ > exp(4.9 — c ) « 80 even in the limit 
A — > oo, a result relevant for 4 He and marginally relevant 
for the high-T c superconductors with k ss 50 - 100. A 
(more physical) alternative to stabilize the twisted-pair 
configuration is to decrease the surrounding pressure via 
reducing the vortex density. For fields B < H Cl the in- 
teraction becomes short-range and we have to return to a 
finite screening length A < oo in V. Indeed, in the dilute 
limit we find a metastable TP for fields B < 1.65 H Cl , i. 
e., a > 1.8 A, where k — 10 has been chosen. 

Next we discuss the twisted-triplet state in the vortex- 
lattice (results for a /£ = 10 are quoted in the text; 
see Table I for a summary). Following the scheme de- 
scribed above we find a metastable TT configuration 
well in the London regime, i. e., with a minimal sep- 
aration drain ~ 0.8 a £ between the vortices. Due 
to the hexagonal symmetry of the surrounding lattice, 
the twist is constrained within a length L TT sa 2.2 a 
and its energy is high as compared to the rectilinear 
groundstate, E TT — E w 2.1e a , see Fig. 2. On the 
other hand, the barrier stabilizing the metastable state 
is small, E s — E TT w 0.32 e a . The saddle still is in 
the London regime with d m i n > 0.3 a . The reswitching 
goes through a first collapse, leading to the creation of a 
small transverse loop which subsquently shrinks to zero, 
leaving behind three rectilinear vortices in their ground- 
state, see the inset of Fig. 2. Note the sharp drop in 
energy when crossing the saddle, a consequence of the 
sudden reswitching and an indication that the present 
constraint does not produce a force-free configuration at 
the saddle-point. We defer the detailed discussion of the 
various types of (re)switching saddles (loop creation, hys- 
tcrctic effects) to a forthcoming paper. 

Model vortex-liquid: Here we reduce the description 
to those vortices directly involved in the entanglement, 
i. e., n = 2 and n = 3 for the twisted-pair and twisted- 
triplet, respectively. We choose the surrounding potential 
mimicking the vortex-liquid to be of the form 

V(R) = {(* 2 + .sV), 

with / = 2/al (/ = 3/a 2 ,) for the TP (TT) configuration; 
s denotes an asymmetry parameter. In their analysis of 
the TP problem, Carraro and Fisher |l3| chose a circular 
symmetry with s = 1, where the metastable twisted state 
takes a helical shape (with a pitch L p ) and is degenerate 
with the rectilinear groundstate in the limit L p — > oo. 



In this limit, the reswitching barrier is characterized by 
a high degree of symmetry and using a /£ = 10 we find 
E a — Eh w 1.0e o a o , in agreement with previous work 
p3[ . However, in a liquid phase one does not expect 
the potential acting on a vortex pair to have a circular 
symmetry; rather the pressure of the surrounding vor- 
tices produces an asymmetric potential with s > 1 (one 
expects s w 2.7 on simple geometrical grounds, assum- 
ing a similar pressure to act in the liquid phase as in 
the solid one). Increasing s beyond the critical value 
s c w 1.3 the entangled configuration becomes unstable 
and no metastable TP configuration is expected to ex- 
ist in the liquid phase either. One could argue, that in 
a vortex liquid the asymmetric potential can accommo- 
date itself and rotate along with the pair, however, such 
a configuration involves the twist of four lines at least. 
The approximation of a circular effective potential mim- 
icking the pressure of the surrounding liquid is already 
quite reasonable for the triplet (and further improves for 
the configurations involving larger loops with 6 and more 
vortices). In fact, the main reason for studying the triplet 
as the elementary entangled configuration is the existence 
of a metastable twisted configuration in the solid, im- 
plying the existence of a similar metastable state in the 
liquid which still exhibits lattice order on short spatial 
and temporal scales; in a (viscous) liquid, the entangle- 
ment will first go through a confined high-energy twist, 
similar to the one in the lattice, which subsequently re- 
laxes into a low-energy helical state due to the absence 
of shear forces in the liquid. The excitation energy of 
this relaxed state depends on the pitch L p of the he- 
lix and we have verified that for L p 3> a the tilt energy 
E h -E fts (27r 2 /9)[ln(a o /£)+0.5]£ o ao/-kp provides a good 
approximation within the London regime (note that here 
the vortices twist by the angle 27r/3 on the length L p ). 
The saddle-point energy for (re)switching remains high, 
however; in the limit L p — > oo an analysis based on sym- 
metry arguments similar to those used by Carraro and 
Fisher |13[ provides the result E s — Eh w 2Ae a , where 
a /£_ = 10 and A — > oo has been chosen. Results for other 
values of a /£ are summarized in Table I. 

In conclusion, we have investigated the elementary en- 
tangled configurations, the twisted-pair and the twisted- 
triplet in a vortex-lattice and -liquid phase. Away from 
the dilute limit, we find that the twisted-pair is unsta- 
ble and hence irrelevant for the discussion of entangle- 
ment. The basic loop of entanglement is the twisted- 
triplet which is metastable both in the vortex-lattice and 
-liquid phase. In the vortex-lattice, the TT is a high en- 
ergy state stabilized by a comparatively small barrier and 
hence the vortex-lattice does not entangle. The high en- 
ergy of the TT excitation is a consequence of the confine- 
ment produced by the lattice symmetry of the surround- 
ing potential. Upon melting, the shear modulus disap- 
pears and the hexagonal lattice symmetry changes to a 
rotational one in the liquid. The high-energy twisted- 



3 



triplet dissolves into a low-energy helix with a pitch L p 
determined by the mutual entanglement in the liquid. 
It is the phenomenon of deconfmement and its associ- 
ated drop in excitation energy which tends to bind the 
two (originally unrelated) transitions of melting (loss of 
translational lattice symmetry) and entanglement (loss of 
longitudinal superconductivity) together. Whereas melt- 
ing immediately triggers entanglement in a system with 
short range interactions, it remains to be shown whether 
a disentangled liquid state can be stabilized in a system 
characterized by an interaction with long range. 
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FIG. 1. Twisted-triplet configuration embedded within a 
vortex-lattice. The 9 nearest neighbors are allowed to re- 
lax when twisting the inner vortex triplet. A confined 
(L T t ~ 2.2 a ), high energy metastable state is found, with an 
excitation energy E TT — E a w 2.12 e a and stabilized against 
reswitching by the small barrier E a ~ 0.32 e a . 



FIG. 2. Excitation energy E T T(d ) — E versus distance 
d for the clamped twisted-triplet configuration in a vor- 
tex-lattice. Inset shows top-view of the clamped configura- 
tions marked by the arrows, particularly the metastable state 
and the fast reswitching geometry close to the saddle. Note 
the small relaxation amplitude of the surrounding vortices. 
Upon melting, the confined high-energy twisted configuration 
turns into a low-energy helical state, whereas the saddle-point 
configuration for (re)switching stays confined and remains at 
high energy (see Table I). 



TABLE I. Numerical results obtained for the triplet con- 
figurations in the vortex-solid (s) and -liquid (1) phases. 
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